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ABSTRACT 


The  Building  and  Fire  Research  Laboratory  (BFRL)  was  given  the  opportunity  to  make  measurements 
during  fire  calibration  tests  of  the  heat  detection  system  in  an  aircraft  hangar  with  a nominal  30.4  m 
(100  ft)  ceiling  height  near  Dallas,  TX.  Fire  gas  temperatures  resulting  from  an  approximately 
8250  kW  isopropyl  alcohol  pool  fire  were  measured  above  the  fire  and  along  the  ceiling.  The  results 
of  the  experiments  were  then  compared  to  predictions  from  the  computer  fire  models  DETACT-QS, 
FPETOOL  and  LAVENT.  In  section  A of  the  analysis  conducted,  DETACT-QS  and  FPETOOL 
significantly  underpredicted  the  gas  temperatures.  LAVENT  at  the  position  below  the  ceiling 
corresponding  to  maximum  temperature  and  velocity  provided  better  agreement  with  the  data.  For 
large  spaces,  hot  gas  transport  time  and  an  improved  fire  plume  dynamics  model  should  be 
incorporated  into  the  computer  fire  model  activation  routines.  A computational  fluid  dynamics  (CFD) 
model,  HARWELL  FLOW3D,  was  then  used  to  model  the  hot  gas  movement  in  the  space. 

Reasonable  agreement  was  found  between  the  temperatures  predicted  from  the  CFD  calculations  and 
the  temperatures  measured  in  the  aircraft  hangar.  In  section  B,  an  existing  NASA  high  bay  space  was 
modelled  using  the  computational  fluid  dynamics  model.  The  NASA  space  was  a clean  room,  27.4  m 
(90  ft)  high  with  forced  horizontal  laminar  flow.  The  purpose  of  this  analysis  is  to  determine  how  the 
existing  fire  detection  devices  would  respond  to  various  size  fires  in  the  space.  The  analysis  was 
conducted  for  32  MW,  400  kW,  and  40  kW  fires. 


Key  words:  clean  rooms;  computer  models;  computational  fluid  dynamics;  detector  response;  field 
modeling;  fire  detection;  fire  models;  forced  air  flow;  fire  plumes;  fire  tests;  high  bays;  response 
time;  sprinkler  response 
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1.0  INTRODUCTION 


Large  spaces,  such  as  those  found  in  warehouses,  historical  buildings,  atriums,  and  aircraft  hangars, 
represent  some  of  the  most  difficult  fire  protection  challenges  since  they  are  frequently  of  historical 
significance,  contain  large  quantities  of  fuel,  have  unique  geometries,  and/or  present  special  life  safety 
problems.  Accurate  detector  activation  predictions  are  important  in  these  large  spaces,  as  timely 
detection  of  a fire  is  more  difficult  due  to  the  distance  heat  and  products  of  combustion  must  travel  to 
reach  sprinklers  and  detectors.  Since  fires  frequently  grow  at  an  exponential  rate,  an  increased  time 
to  detection  results  in  larger  fires  at  the  time  of  detection  and  larger  fires  to  be  suppressed  by,  e.g., 
an  automatic  sprinkler  system.  Even  a modest  uncertainty  in  the  prediction  of  the  activation  time  may 
lead  to  a large  uncertainty  in  the  fire  size  used  to  predict  the  hazard,  a central  element  of  many  fire 
protection  analyses. 

There  has  been  substantial  effort  to  verify  activation  predictions  in  small  and  medium  sized  rooms, 
but  little  in  large  spaces.  In  January,  1993,  a report  describing  model  verification  tests  conducted  in  a 
14  m (46  ft)  high  hangar  was  released  by  NIST  [1].  Although  the  correlations  in  the  computer  fire 
models  were  based  on  ceiling  heights  up  to  16  m (52  ft),  previous  to  this  study,  comparisons  of 
computer  fire  model  predictions  with  experimental  measurements  for  ceiling  heights  approaching  that 
height  had  not  been  made.  Conducting  verification  experiments  for  large  spaces  is  difficult  due  to  the 
lack  of  availability  of  adequate  facilities  for  live  fire  tests.  The  largest  dedicated  fire  test  facility  in 
the  United  States  is  limited  to  a height  of  18.2  m (60  ft). 

Why  is  model  verification  important?  It  is  the  responsibility  of  the  user  of  a model  to  determine  the 
suitability  of  a particular  fire  model  for  a given  situation.  A successful  comparison  of  the  model 
predictions  with  real-scale  fire  results  helps  the  user  make  that  decision.  Comparison  of  a fire  model 
against  appropriately  chosen  experiments  can  provide  the  necessary  base  to  extend  the  use  of  the 
model  to  real  fire  applications.  Comparisons  also  allow  model  developers  to  examine  the  accuracy  of 
a model  and  identify  the  needs  for  future  model  development.  It  is  critical  that  the  user  have 
comparisons  of  model  predictions  to  a wide  range  of  experimental  data. 

The  Building  and  Fire  Research  Laboratory  (BFRL)  was  given  the  opportunity  to  make  measurements 
during  fire  calibration  tests  of  the  heat  detection  system  in  an  aircraft  hangar  with  a nominal  30.4  m 
(100  ft)  ceiling  height  near  Dallas,  TX.  Three  closed-door  tests  were  conducted. 

This  project  was  funded  in  two  phases.  The  U.S.  General  Services  Administration  (GSA)  provided 
initial  funding  to  conduct  the  experiments  and  to  compare  measured  data  with  the  predictions  of 
DETACT-QS,  LAVENT,  and  FPETOOL  computer  fire  models.  GSA’s  interest  is  to  verify  and 
provide  recommendations  for  the  use  of  the  GSA  Engineering  Fire  Assessment  System  for  predicting 
the  response  of  sprinklers  and  detectors  in  large  spaces.  FPETOOL  is  the  main  component  of  this 
system.  The  National  Aeronautics  and  Space  Administration  (NASA)  provided  funding  to  compare 
the  predictions  of  a computational  fluid  dynamics,  or  field  model  to  the  experiment.  Comparisons  of 
the  experimental  data  to  both  the  computer  fire  models  and  the  field  model  are  described  in  section  A 
of  this  report.  NASA  also  funded  section  B of  this  project  which  uses  the  field  model  to  predict 
smoke  movement  in  a NASA  Goddard  high  bay  space.  NASA’s  interest  is  in  the  ability  of  field 
modelling  to  provide  information  that  would  be  useful  in  selecting  effective  fire  protection  options  for 
their  high  bay  spaces. 
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2.0  EXPERIMENTAL 


2.1  Building  layout 

The  aircraft  hangar,  pictured  in  figure  1,  measured  approximately  389  m (1280  ft)  long  by  115  m 
(377  ft)  deep  (81  m (270  ft)  to  a firewall),  and  30.4  m (100  ft)  high.  TTie  hangar  had  seven  bays  and 
was  capable  of  housing  seven  wide-bodied  aircraft  side  by  side.  The  hangar  contained  draft  curtains 
which  were  spaced  approximately  12.5  m (41  ft)  apart  as  shown  in  figure  2.  Draft  curtains  sue 
partitions  extending  vertically  down  from  and  perpendicular  to  the  ceiling.  Draft  curtains  are 
designed  to  slow  the  spread  of  combustion  gases  across  the  ceiling.  The  hot  gases  collect  between  the 
draft  curtains  and  flow  downward  until  spilling  into  the  adjacent  curtained  areas.  This  concentration 
of  smoke  and  hot  gases  inside  the  draft  curtain  is  intended  to  speed  up  the  operation  of  sprinklers 
within  the  curtained  area.  Draft  curtains  in  this  facility  were  3.7  m (12  ft)  deep.  A plan  view  of  the 
fire  test  curtained  area  and  the  two  adjacent  curtained  areas  is  shown  in  figure  2. 

The  ceiling  geometry  is  shown  in  figure  3.  Perpendicular  to  the  draft  curtains,  a series  of  large  I- 
beams  are  located  near  the  ceiling  which,  when  coupled  with  the  draft  curtain  geometry,  appear  to 
partition  the  curtained  areas  into  cubes.  A second  set  of  I-beams  runs  both  below  the  draft  curtains 
and  the  I-beams  forming  a second  set  of  open  cubes  at  a height  of  24  m (79  ft)  off  the  floor.  Large  I- 
beams  join  the  top  and  bottom  I-beams  in  a triangular  structure.  The  underside  of  the  roof  is  made  of 
corrugated  steel  with  38  mm  (1.5  in)  deep  corrugations. 

The  front  of  the  hangar  consisted  primarily  of  sliding  doors,  which  if  ail  were  open  would  provide  a 
maximum  opening  289  m (948  ft)  in  length  and  22  m (72  ft)  high.  The  hangar  also  contained  many 
personnel  doors  both  to  the  outside  and  to  other  parts  of  the  building.  The  large  doors  were  closed 
during  the  tests;  however,  some  of  these  personnel  doors  may  well  have  been  open  during  the  tests. 
(There  was  much  work  going  on  in  the  building  during  testing.)  The  building  HVAC  system  was  off 
for  the  duration  of  the  three  tests. 

2.2  Test  fire 

The  test  fire  was  an  array  of  nine  pans,  each  0.91  m x 0.91  m (3  ft  x 3 ft),  placed  together  to  form  a 
square  fire  area  totalling  7.5  m^  (81  ft^)  as  shown  in  figure  4.  The  effective  circular  diameter  of  this 
fire  is  3.1  m (10  ft).  The  pans  were  0.076  m (0.25  ft)  deep.  The  pans  were  supported  on  bricks, 

127  mm  (0.42  ft)  off  the  floor.  The  test  fuel  was  technical  grade  isopropyl  alcohol.  A total  of  135 
liters  (35.1  gal)  of  fuel  was  used  per  test,  15  liters  (4.0  gal)  of  fiiel  per  pan.  The  fiiel  depth  in  each 
pan  was  approximately  13  mm  (0.5  in).  Cloth  wicks  were  draped  between  the  pans  to  aid  in  fire 
spread  after  ignition  of  one  pan. 

The  fire  was  located  on  the  floor  in  the  center  of  the  hangar  building.  At  8 seconds  after  ignition  of 
the  first  pan,  all  the  fiiel  pans  were  fiilly  involved.  A picture  of  the  fire  at  steady-state  is  shown  in 
figure  5.  Steady  burning  was  maintained  for  approximately  three  and  one-half  minutes,  after  which 
the  flame  height  started  to  decrease.  The  fire  was  allowed  to  bum  out,  which  took  several  additional 
minutes. 

As  there  was  no  load  cell  or  oxygen  calorimeter  available  for  use  during  these  experiments,  the  heat 
release  rate  of  the  fire  was  estimated  from  previous  bum  tests  conducted  at  NIST.  These  tests  are 
described  in  [1].  Isopropyl  alcohol  pan  fires  3.3  m^  (36  ft^)  were  conducted  under  the  NIST  large 
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calorimeter.  These  tests  used  4-0.91  x 0.91  m (3  x 3 ft)  steel  pans.  The  effective  circular  diameter 
of  this  fire  was  2.1  m (6.9  ft).  In  the  bum  tests  conducted,  the  average  burning  rate  was  0.036 
kg/m“/s  (0.0074  Ib/fr/s),  and  the  average  heat  release  rate  was  3650  kW  (3460  Btu/s).  The  net  heat 
of  combustion  is  reported  in  the  NFPA  handbook  [2]  as  30.45  MJ/kg  (13107  Btu/lb). 

Assuming  there  is  sufficient  oxygen  available  for  complete  combustion,  the  heat  release  rate  of  a fire 
can  be  calculated  as: 

Q = lii"  * Ahj.  * Af 

where: 

Q = heat  release  rate  (kW) 
m"  = burning  rate  (kg/m^/s) 

Ahj.  = net  heat  of  combustion  (kJ/kg) 

Af  = surface  area  of  fuel  (m-) 

To  characterize  the  burning  rate  of  a liquid  pool  fire,  it  is  required  that  a classification  be  made 
according  to  the  dominant  heat  transfer  mechanism,  which  can  be  dependent  on  pool  diameter  [3]. 
When  the  pool  diameter,  D,  is  very  small,  D < 0.03  m (0.01  ft),  conductive  heat  transfer  determines 
the  rate  of  burning,  while  the  radiative  term  predominates  if  D is  large,  D > 1.0  m (3.3  ft).  At  D = 
1.0  m (3.3  ft)  and  above,  the  burning  rate  has  been  shown  to  be  constant  [4].  Since  the  isopropyl 
alcohol  pan  fires  conducted  at  NIST  were  for  an  effective  diameter  of  2.1  m (6.9  ft),  which  is  well 
within  the  radiative  dominated  range,  the  same  burning  rate  was  used  for  the  3.1  m (10  ft)  effective 
circular  diameter  fire  conducted  in  the  aircraft  hangar.  Using  the  above  equation,  the  calculated  heat 
release  rate  for  the  fire  in  the  aircraft  hangar  is  approximately  8250  kW  (7820  Btu/s). 

The  heat  release  rate  as  a function  of  time  is  shown  in  table  1.  An  initial  growth  phase  of  the  fire 
was  selected  to  provide  a reasonable  fit  with  the  measured  temperature  rise  early  in  the  fire,  as  well 
as  with  test  observations. 

The  fraction  of  the  heat  release  rate  that  is  radiated  from  the  flaming  region  directly  affects  the  gas 
temperature.  The  radiant  heat  fraction  is  sometimes  referred  to  as  X^.  For  many  common  fuels  with 
luminous  flames,  it  has  been  found  that  the  radiative  fraction  is  35  percent  [5].  Although  many 
alcohols  bum  without  luminous  flames,  isopropyl  alcohol  flames  are  very  luminous,  appearing  much 
like  the  flames  produced  from  burning  heptane.  Since  isopropyl  alcohol  bums  with  luminous  flames, 
the  analysis  was  performed  with  a radiative  fraction  of  35  percent. 

To  determine  the  sensitivity  of  the  gas  temperature  predictions  to  heat  release  rate,  a 20  percent 
higher  heat  release  rate,  9900  kW  (9375  Btu/s)  was  also  examined  in  this  report.  This  provides  a 
higher  end  prediction  of  the  temperatures  in  the  hangar.  The  higher  heat  release  rate  analysis  takes 
into  account  possible  uncertainties  in  measurement  of  prior  tests  of  the  burning  rate  performed  in  the 
NIST  large  calorimeter. 
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2.3  Measurements 


2.3.1  Type  of  measurements  taken 

Thermocouples  were  used  to  measure  the  gas  temperatures  above  the  fire  and  along  the  ceiling.  The 
thermocouples  were  bare-bead,  0.52  mm  (0.02  in)  diameter,  type  K,  chromel-alumel  thermocouples. 
Brass  disks  were  used  to  simulate  the  heating  of  fusible  links.  The  disks  were  9.5  mm  (0.38  in) 
diameter,  0.81  mm  (0.03  in)  thick  brass.  The  RTI  of  the  brass  disks  is  33.4(m-s)^'^^  (60.5  (ft-s)^^^). 
Figure  6 shows  the  instrumentation  layout.  The  temperatures  of  the  thermocouples  and  brass  disks 
were  recorded  every  4 seconds  on  a computer-based  data  acquisition  system. 

To  measure  the  gas  temperature,  six  temperature  measurements  were  taken  vertically  below  the 
ceiling  east  of  the  fire  at  radial  distances  of  1.5,  3.0,  4.6,  and  6.1  m (5,  10,  15,  and  20  ft)  from  the 
centerline  of  the  fire  and  west  of  the  fire  at  a radial  distance  of  3.0  m from  the  centerline  of  the  fire. 
At  each  radial  position,  a thermocouple  and  calibrated  brass  disk  were  placed  at  0. 15  m (0.5  ft)  and 
0.91  m (3  ft)  below  the  ceiling.  Thermocouples  were  also  placed  at  0.30  m (1.0  ft)  and  1.5  m (5  ft) 
below  the  ceiling.  Assuming  that  the  smoke  and  temperature  flow  together,  the  thermocouple 
readings  can  indicate  the  initial  arrival  of  the  smoke  to  the  thermocouple  location.  This  assumption  is 
supported  by  the  experiments  of  Heskestad  and  Delichatsios  [6]. 

To  track  smoke  filling  in  the  fire  test  curtained  area,  thermocouples  were  placed  to  the  east  and  west 
of  the  fire  at  5.5  m (18  ft)  radially  out  from  the  fire  centerline.  Thermocouples  were  placed 
vertically  from  the  ceiling  to  the  bottom  of  the  draft  curtain  as  shown  in  figure  6 to  determine  the 
layer  height. 

To  measure  the  centerline  plume  temperatures,  eight  thermocouples  were  placed  directly  over  the  fire 
at  0.025,  3.0,  6.1,  9.1,  12.2,  15.2,  18.2,  and  21.3  m (0.082,  10,  20,  30,  40,  50,  60,  and  70  ft) 
below  the  ceiling. 

To  track  the  flow  of  smoke,  a thermocouple  was  placed  0.025  m (0.082  ft)  directly  under  the  draft 
curtains  between  the  test  area  and  the  adjacent  curtained  areas.  An  array  of  four  thermocouples  was 
placed  at  the  center  of  the  east  adjacent  curtained  area  and  of  the  west  adjacent  curtained  area. 
Thermocouples  were  also  placed  0.025  m (0.082  ft)  below  the  draft  curtains  at  the  far  end  of  the  east 
adjacent  curtained  area  and  at  the  far  end  of  the  west  adjacent  curtained  area.  Thermocouples  were 
placed  to  the  north  and  south  of  the  fire,  at  radial  distances  of  3.0,  6.1,  and  9.1  m from  the  pan 
centerline  (10,  20,  and  30  ft),  each  0.15  m (0.5  ft)  below  die  ceiling, 

2.3.2  Uncertainty  of  measurements 

The  manufacturer  of  the  thermocouple  wire  [7]  states  that  there  is  a ± 2.2  °C  uncertainty  in  the 
thermocouple  reading.  The  resolution  of  the  data  acquisition  system  is  0.3  °C.  Experimental 
uncertainties  such  as  re-radiation  and  chemical  dissociation  are  negligible  for  the  temperature  range 
realized  in  this  experiment. 

The  fire  pulsations  are  reflected  in  the  temperature  fluctuations  seen  in  the  curves  for  measured 
temperature.  A mean  temperature  for  each  measurement  point  was  calculated  and  is  reported  with  the 
associated  standard  deviation.  This  deviation  is  plus  or  minus  three  to  seven  degrees  Celsius 
depending  on  the  location  of  the  thermocouple. 
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3.0  COMPUTER  MODELS 


The  NFPA  handbook  groups  computer  fire  models  into  two  categories;  enclosure  fire  models  and 
special-purpose  models  [8].  Enclosure  fire  models  are  structured  to  predict  the  interaction  of  multiple 
fire  processes  involving  heat  transfer,  fluid  mechanics,  and  combustion  chemistry  occurring  in  an 
enclosure.  Special  purpose  computer  models  are  computer  models  designed  for  special  purpose 
analyses  such  as:  structural  fire  resistance,  prediction  of  response  time  of  heat  detectors  and  automatic 
sprinklers,  egress  design,  or  performance  of  smoke  control  or  ventilation  systems.  Fire  models  are 
further  subdivided  into  probabilistic  and  deterministic  models.  Probabilistic  models  treat  fire  growth 
as  a series  of  sequential  events  or  states.  Probabilities  are  assigned  to  each  transfer  point  based  on 
analysis  of  relevant  experimental  data  and  historical  fire  incident  data.  In  contrast,  deterministic 
models  represent  the  processes  encountered  in  a compartment  fire  by  interrelated  mathematical 
expressions  based  on  physics  and  chemistry.  According  to  this  system  classification,  DETACT-QS, 

LA  VENT,  and  the  activation  subroutine  in  FPETOOL  are  all  special  purpose  deterministic  fire 
models  [8]. 

A CFD  model  can  be  used  as  a deterministic  enclosure  model.  A CFD  model  solves  the  mass, 
energy,  and  momentum  equations  directly  which  provides  more  flexibility  in  modelling  complex 
situations.  In  using  a CFD  model,  the  enclosure  is  divided  up  into  many  small  spaces  or  grid 
volumes  and  the  conservation  equations  are  then  solved  to  yield  the  values  of  temperature,  velocity, 
etc.  at  the  center  of  these  grid  volumes.  CFD  modelling  then  provides  temperature  and  velocity 
distributions  across  the  enclosure. 

3.1  DETACT-QS  and  FPETOOL 

It  is  important  to  understand  the  physics  of  the  computer  models  and  the  assumptions  built  into  each 
code.  FPETOOL,  written  by  Nelson  [9,10]  and  DETACT-QS  (DETector  ACTuation  Quasi-Steady), 
written  by  Evans  and  Stroup  [11,12],  are  each  based  on  experimental  correlations  developed  by 
Alpert  for  steady-state  fires  [13].  These  correlations  give  the  maximum  temperature  and  maximum 
velocity  of  the  ceiling  jet  as  a function  of  the  steady  state  heat  release  rate  of  the  fire,  the  radial 
distance  from  the  fire  to  the  detector,  and  the  height  of  the  ceiling  above  the  fuel.  Both  programs 
assume  the  detector  or  sprinkler  is  located  at  a distance  below  the  ceiling  corresponding  to  maximum 
temperature  and  velocity.  Alpert’s  correlations  assume  a smooth,  unconfined  ceiling.  They  also 
assume  that  steady-state  correlations  can  be  applied  to  a growing  fire  over  small  time  intervals.  Both 
DETACT-QS  and  FPETOOL  perform  a linear  interpolation  to  determine  the  heat  release  rate  at 
intermediate  times. 

In  both  programs,  the  transport  time  of  the  smoke  and  hot  gases  from  the  fire  to  the  thermal  detector 
is  neglected.  Also  in  both  programs,  the  detector  is  subject  to  the  maximum  temperature  and  velocity 
of  the  ceiling  Jet.  FPETOOL  accounts  for  the  impact  of  the  hot  gases  entrained  into  the  ceiling  Jet  on 
the  temperature  and  velocity  of  the  Jet  as  it  passes  through  the  hot  smoke  layer;  DETACT-QS  does 
not.  The  expressions  used  in  DETACT-QS  and  FPETOOL  for  temperature  and  velocity  are  based  on 
a nondimensional  parameter  r/h,  where  r is  the  radial  distance  from  the  centerline  of  the  fire  to  the 
thermal  device,  and  h is  the  height  from  the  base  of  the  fire  to  the  ceiling.  These  expressions  are 
independent  of  radius  for  r/h  < 0.18,  which  is  assumed  to  be  within  the  plume  region. 

The  user  inputs  for  DETACT-QS  are  shown  in  table  2.  The  user  inputs  for  FPETOOL  are  shown  in 
table  3.  FPETOOL  version  3.01  was  used.  In  both  programs,  999  °C  is  used  for  the  response 
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temperature  of  the  device  as  the  predicted  thermal  device  temperature  was  compared  to  the  measured 
temperature  of  the  disk  and  not  to  the  activation  time  of  a sprinkler  or  detector,  i.e.,  the  device  was 
programmed  not  to  go  off.  FPETOOL  and  DETACT-QS  do  not  model  draft  curtains.  DETACT-QS 
has  no  walls  and  assumes  open  burning.  It  can  not  accommodate  heat  losses  to  the  ceiling.  To  model 
a draft  curtain  using  FPETOOL,  the  dimensions  of  the  curtained  area  are  used  for  the  room  spacial 
dimensions,  the  program  assumes  there  are  walls.  The  door/window  option  is  then  used  to  simulate 
the  area  under  the  draft  curtains  as  large  openings  in  the  walls. 

3.2  LAVENT 

LAVENT  (Link-Actuated  VENT),  is  a program  written  by  Davis  and  Cooper  [14,15],  to  simulate  the 
environment  and  response  of  sprinkler  links  in  compartment  fires  with  draft  curtains  and  fusible-link- 
actuated  ceiling  vents.  LAVENT  is  similar  to  FPETOOL  and  DETACT-QS  in  that  it  assumes  that 
steady  state  correlations  can  be  applied  to  a growing  fire  over  small  time  intervals;  and  it  too  neglects 
the  transport  time  of  the  smoke  and  hot  gases  from  fire  to  thermal  detector.  LAVENT  does  account 
for  the  impact  of  the  hot  upper  layer  on  the  ceiling  jet.  The  important  difference  between  LAVENT 
and  the  other  activation  models,  is  that  LAVENT  accounts  for  position  of  the  thermal  detector  below 
the  ceiling  in  the  ceiling  jet.  It  cannot  model  recessed  or  concealed  sprinklers,  as  the  velocity  in 
LAVENT  is  zero  at  the  ceiling.  The  expressions  for  ceiling  jet  temperature  and  velocity  used  in 
LAVENT  are  independent  of  radius  for  r/h  < 0.2,  which  is  assumed  to  be  within  the  plume  region. 
The  expression  for  ceiling  jet  temperature  and  velocity  as  a ftmction  of  height  beneath  the  ceiling  used 
in  LAVENT  are  not  valid  in  the  plume  region,  where  there  are  no  appreciable  temperature  or  velocity 
gradients  below  the  ceiling.  LAVENT  will  only  provide  an  accurate  temperature  prediction  in  the  r/h 
0.2  region  at  the  point  below  the  ceiling  where  the  maximum  temperature  and  velocity  occur.  This 
unfortunately  must  be  determined  by  the  user  through  trial  and  error,  although  some  guidance  is 
provided  in  the  user  guide. 

LAVENT  utilizes  draft  curtains  instead  of  walls.  LAVENT  models  just  one  curtained  area.  When 
that  area  fills  with  smoke  it  spills  out  into  an  infinite  open  space.  The  user  inputs  for  LAVENT  are 
shown  in  table  4. 

3.3  FLOW3D 

A computational  fluid  dynamic  (CFD)  or  field  model,  HARWELL  FLOW3D  [16],  was  used  to 
simulate  the  fires  in  the  large  spaces  analyzed  in  this  report.  This  model  uses  a finite  difference 
method  to  solve  the  Navier  Stokes  equations  in  either  two  or  three  dimensions  on  a user-specified 
grid.  The  Navier-Stokes  equations  include  mass,  momentum,  and  energy  conservation. 

While  the  CFD  model  enjoys  a degree  of  flexibility  that  is  missing  in  the  special  purpose  activation 
models,  there  are  a few  drawbacks  to  using  a CFD  model.  The  model  is  typically  more  difficult  to 
use  in  that  an  appropriate  grid  must  be  chosen  to  represent  the  physical  space,  accurately  solve  the 
partial  differential  equations,  and  provide  for  adequate  representation  at  the  boundaries.  These 
problems  typically  require  thousands  of  grid  volumes  to  provide  accurate  results.  Long  computation 
times  are  therefore  required  for  CFD  calculations.  Fire-based  CFD  calculations  require  a turbulence 
model  in  order  to  simulate  the  gas  momentum  and  energy  flow  in  turbulent  regions.  The  turbulence 
models  now  available  either  do  not  do  an  adequate  job  modelling  fire  plumes  or  have  not  been 
thoroughly  tested  over  a large  range  of  fires. 
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For  the  calculations  reported  here,  the  fire  is  represented  by  a heat  release  rate  with  the  energy  input 
being  distributed  over  several  grid  volumes  to  avoid  an  unrealistically  large  energy  input  in  any  single 
grid  volume.  The  ideal  gas  law  is  used  to  represent  the  equation  of  state  for  air.  Buoyancy  is 
incorporated  in  the  calculation  by  assuming  that  the  gas  is  compressible  and  that  the  density  is  given 
by  the  ideal  gas  law.  Fire  induced  turbulence  is  modeled  using  the  k-e  model  16].  The  turbulence 
model  depends  on  a number  of  user  specified  parameters  which  are  given  in  table  5.  Wall  boundary 
conditions  are  treated  using  the  standard  model  for  turbulent  flow.  Table  6 lists  the  key  words  used 
in  the  command  file  to  activate  physical  models,  the  grid  volume  dimensions  for  some  of  the 
calculations,  the  ambient  variables,  and  the  key  words  to  specify  the  boundary  conditions  used  for  the 
FLOW3D  calculations. 


4.0  COMPARISON  OF  COMPUTER  FIRE  MODELS  TO  MEASUREMENTS 

4. 1 Plume  symmetry 

Table  7 shows  the  measured  mean  temperature  and  standard  deviation  of  the  gas  at  3.0  m (10  ft)  to 
the  east  of  the  pan  centerline  and  3.0  m (10  ft)  to  the  west  of  the  pan  centerline  at  a depth  of  0.15  m, 
0.30  m,  0.91  m,  and  1.5  m (0.5  ft,  1 ft,  3 ft,  and  5 ft)  below  the  ceiling  for  the  three  experiments 
conducted.  By  comparing  the  east-west  temperature  measurements,  it  can  be  seen  that  the  smoke 
plume  was  most  symmetrical  in  test  2.  Based  on  plume  symmetry,  test  2 was  used  for  comparison  to 
the  predictions  of  the  computer  models. 

4.2  Gas  temperatures 

Table  8 shows  the  steady  state  gas  temperatures  at  0,  1.5,  3.0,  4.6,  6.1,  and  12.5  m (0,  5,  10,  15, 

20,  and  41  ft)  radial  distance  from  the  centerline  of  the  fire  as  measured  in  the  experimental  section. 
The  measured  temperature  reported  is  the  mean  temperature  from  32  to  100  s.  The  number  in 
parentheses  following  the  mean  temperature  is  the  standard  deviation.  This  time  period  represents 
quasi-steady  conditions  of  the  temperature  near  the  ceiling. 

4.2.1  Inside  the  fire  test  curtained  area 

Figures  7 and  8 show  the  graphical  comparison  of  the  measured  and  predicted  gas  temperatures  at  a 
radial  distance  of  1 .5  m (5  ft)  and  6. 1 m (20  ft)  from  the  centerline  of  the  fire  respectively.  These 
radial  distances  correspond  to  non-dimensional  radii  (r/h)  of  0.05,  and  0.2.  The  measured  gas 
temperatures  are  shown  as  solid  lines.  In  each  graph,  these  four  solid  lines  represent  the  measured 
temperatures  at  0.15  m,  0.30  m,  0.91  m,  and  1.5  m,  (0.5  ft,  1 ft,  3 ft,  and  5 ft)  below  the  ceiling. 
Since  the  temperatures  at  these  four  distances  below  the  ceiling  are  coincident,  they  were  not 
individually  labeled.  These  curves  fluctuate,  which  reflects  the  pulsation  of  the  fire.  The  predictions 
of  DETACT-QS  and  FPETOOL  are  also  shown  as  solid  lines.  The  predicted  values  are  smooth  and 
straight.  Predictions  of  gas  temperature  by  LAVENT  are  shown  as  various  style  dashed  lines. 

Arrows  on  the  graph  point  to  predictions  of  LAVENT  at  vertical  distances  below  the  ceiling  of  d = 
0.15  m,  0.30  m,  0.91  m,  and  1.5  m (0.5  ft,  1 ft,  3 ft,  and  5 ft).  There  is  only  one  curve  shown  for 
DETACT-QS  and  FPETOOL  at  each  radial  distance  as  these  programs  do  not  account  for  the  vertical 
distance  below  the  ceiling. 
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Measurement  locations  from  the  fire  centerline  through  6. 1 m (20  ft)  radial  distance  from  the  fire  are 
within  the  fire  test  curtained  area  as  shown  in  figure  6.  For  these  locations,  the  r/h  values  are  all  0.2 
or  less.  For  r/h  < 0.2,  the  temperature  predictions  of  the  computer  fire  models  do  not  change. 
DETACT-QS  and  FPETOOL  underpredict  the  excess  temperature  (temperature  rise  above  the  ambient 
value  of  22  °C)  by  34  - 41  percent  (depending  on  the  radial  distance  from  the  fire  centerline)  when 
compared  to  the  measured  values  averaged  over  a 32  - 100  second  time  interval.  It  can  also  be 
observed  from  these  graphs,  that  the  predictions  of  LAVENT  nearer  to  the  ceiling  are  closer  to  the 
measured  values,  and  thus  more  accurate  than  the  DETACT-QS  or  FPETOOL  values.  However,  for 
positions  further  from  the  ceiling,  the  predictions  of  LAVENT  are  less  accurate  than  DETACT  or 
FPETOOL.  This  is  due  to  the  measurement  locations  being  inside  the  plume.  The  correlations  for 
ceiling  jet  velocity  and  temperature  as  a function  of  distance  below  the  ceiling  used  in  LAVENT  are 
not  valid  inside  the  plume,  and  thus  the  user  must  be  careful  not  to  use  this  parameter  for  locations 
within  the  plume.  This  phenomena  is  illustrated  in  figure  9. 

However,  for  high  bay  areas,  often  the  location  of  interest  is  an  r/h  < 0.2  due  to  the  ceiling  height. 
Therefore,  when  running  LAVENT  inside  the  plume,  the  user  should  use  only  the  prediction  at  the 
position  below  the  ceiling  corresponding  to  the  maximum  temperature  and  velocity.  This  position 
must  be  found  through  trial  and  error,  although  some  guidance  is  provided  in  [14].  LAVENT  at  the 
point  below  the  ceiling  where  the  temperature  is  maximum  (approximately  0.15  m below  the  ceiling 
in  this  case)  underpredicts  the  excess  temperature  by  6 - 16  percent  when  compared  to  the  measured 
values  averaged  over  a 32  - 100  second  time  interval.  Future  versions  of  LAVENT  should  be 
upgraded  so  that  the  computer  program  does  not  calculate  for  a ceiling  jet  (temperature  gradient  in  the 
vertical  direction)  for  r/h  < 0.2. 

Considering  the  above  mentioned  problems  with  the  predictions  of  the  computer  fire  models,  it  may 
be  useftil  to  use  the  plume  correlations  for  the  plume  centerline  temperatures  to  determine  the 
temperature  rise  above  ambient  for  locations  inside  the  plume  region.  These  plume  correlations 
assume  burning  in  the  open  (i.e.,  no  ceiling),  and  thus  free  entrainment.  There  are  a number  of 
prediction  methods  available  for  plume  temperature,  in  general  producing  very  similar  results.  For 
comparisons  performed  here,  the  prediction  methods  of  McCaffrey  were  selected  [18] 
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where: 

ATg  = centerline  temperature  rise  above  ambient  (K) 

Tg  = ambient  temperature  (K) 
g = acceleration  due  to  gravity  (m/s^) 

C = constant,  0.9 
z = height  above  base  of  fire  (m) 

Q = heat  release  rate  of  the  fire  (kW) 

k,  and  t;  are  constants  whose  value  depends  on  weather  the  location  of  interest  is  within  the  flame, 
intermittent,  or  plume  region  of  the  fire.  This  is  determined  by  the  factor  z/Q^^^  (m/kW^^^).  If  z/Q^^^ 
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is  < 0.08,  it  is  considered  that  this  location  is  in  the  flaming  region  of  the  fire.  In  this  region,  k = 
6.8  m’^'^/s,  and  7;  = 1/2.  If  is  between  0.08  - 0.2,  it  is  considered  that  this  location  is  in  the 
intermittent  region  of  the  fire.  In  the  intermittent  region,  k = 1.9  m/kW‘^^  s,  and  rj  = 0.  If  is 
> 0.2,  it  is  considered  that  this  location  is  in  the  plume  region  of  the  fire  as  is  the  case  for  the 
calculations  near  the  ceiling.  In  the  plume  region  of  the  fire,  k = 1.1  m'^^^/kW^^^.s  and  77  = -1/3. 

The  plume  theory  of  McCaffrey  predicts  a temperature  gradient  from  0.025  m to  1.5  m below  the 
ceiling  of  only  2 °C,  which  compares  well  with  the  measured  temperatures  which  were  coincident 
within  the  uncertainty  band.  The  predicted  temperature  at  0.025  m below  the  ceiling  is  53  °C  and  at 
1.5  m below  the  ceiling  is  55  °C.  The  predicted  temperature  at  0.025  m below  the  ceiling  represents 
an  underprediction  of  the  excess  temperature  by  21  percent  when  compared  to  the  measured  value 
averaged  over  a 32  - 100  second  time  interval. 

Figure  10  shows  the  measured  gas  temperatures  along  the  ceiling,  0.15  m (0.5  ft)  below  the  ceiling  at 
radial  distances  of  0,  1.5,  3.0,  and  4.6  m (0,  5,  10,  and  15  ft).  It  can  be  seen  that  there  is  no 
measured  temperature  gradient  within  the  plume. 

4.2.2  Outside  the  fire  test  curtained  area 

The  measurement  location  at  a radial  distance  of  12.5  m (41  ft)  from  the  centerline  of  the  fire, 
corresponding  to  a nondimensional  radius  of  r/h  = 0.41,  is  outside  the  plume  region  and  on  the  other 
side  of  the  fire  test  curtained  area.  None  of  the  computer  programs  account  for  the  flow  around  the 
draft  curtains,  thus  this  location  can  not  be  modeled  accurately  with  present  computer  fire  models. 

To  model  the  draft  curtain,  and  the  flow  around  the  curtain,  a computational  fluid  dynamics  model  is 
necessary. 

4.2.3  Sensitivity  to  Heat  Release  Rate 

As  was  discussed  in  section  2.2  of  this  report,  an  analysis  was  conducted  to  determine  the  sensitivity 
of  model  predictions  to  the  heat  release  input.  Temperature  predictions  of  the  models  were 
determined  using  a 20  percent  higher  heat  release  rate  of  9900  kW  (9380  Btu/s).  Figure  1 1 shows 
the  DETACT-QS,  FPETOOL,  and  LAVENT  predictions  for  the  higher  heat  release  rate  at  the  6.1  m 
radial  distance  compared  to  the  experimental  data.  This  figure  can  be  compared  to  figure  8 which 
shows  the  temperatures  predicted  for  the  same  location  with  a heat  release  rate  of  8250  kW. 
FPETOOL  and  DETACT-QS  predicted  a 9 percent  (2  °C)  higher  excess  temperature  for  the  20 
percent  higher  heat  release  rate.  LAVENT  predicted  an  11  - 13  percent  (2-4  ° C)  increase  in  excess 
temperature.  The  higher  heat  release  rate  did  change  the  temperature  predictions;  however,  it  did  not 
change  the  overall  results.  That  is,  DETACT-QS  and  FPETOOL  still  underpredicted  the  gas 
temperatures.  LAVENT  predicted  the  temperatures  well  (closer  to  the  experiment)  at  the  position 
below  the  ceiling  corresponding  to  the  maximum  temperature  and  velocity. 

4.3  Hot  gas  transport  time 

The  measured  gas  temperatures  at  the  pan  centerline,  1.5  m,  3.0  m,  4.6  m (5  ft,  10  ft,  15  ft)  and  out 
to  12.5  m (41  ft)  to  the  east  all  show  a significant  jump  in  temperature  at  the  20  s data  point  (table  9 
givestime  to  first  response  of  the  thermocouples)  that  had  not  registered  at  the  16  s data  point.  This 
tells  us  that  the  smoke  plume  reached  these  points  18  ± 2 seconds  after  ignition.  It  is  impossible  to 
determine  more  exactly  than  the  resolution  of  the  data  recorder.  Also,  the  thermocouple  3.0  m (10  ft) 
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to  the  west  of  the  pan  centerline  registered  an  increase  in  temperature  20  seconds  after  ignition.  The 
thermocouple  measurements  combined  with  visual  test  observations,  suggest  the  plume  width  to  be 
approximately  16  m (52  ft).  The  calculated  plume  width  (using  Heskestad’s  plume  theory)  at  the 
ceiling  is  15  m (49  ft)  which  agrees  with  this  result.  [19] 

Figure  12  shows  the  predicted  gas  temperature  by  LAVENT  at  a radial  distance  of  4.6  m (15  ft)  and 
0.15  m (0.5  ft)  below  the  ceiling.  This  location  was  chosen  because  the  measured  temperature 
matched  the  predicted  temperature  so  closely.  Also  shown  in  figure  12  is  the  same  gas  temperature 
prediction  shifted  18  seconds  to  the  right.  The  18  second  shift  was  selected  to  match  the  times  of  the 
predicted  gas  temperature  rise  to  the  times  of  the  measured  temperature  rises.  None  of  the  computer 
fire  models  account  for  the  transport  time  in  the  plume  so  changes  in  the  input  fire  immediately  affect 
the  predicted  temperatures  at  the  ceiling.  A similar  shift  could  be  applied  to  the  FPETOOL  and 
DETACT-QS  results. 


5.0  COMPARISON  OF  FIELD  MODEL  (CFD)  CALCULATIONS  TO  MEASUREMENTS 

The  30.4  m (100  ft)  hangar  curtained  area  was  modelled  using  a heat  source  to  simulate  the  pan  fire. 
Computer  speed  effectively  limits  the  total  number  of  grid  volumes  available  to  model  the  problem. 

It  was  assumed  that  all  the  heat  was  released  in  the  sixteen  model  grid  volumes  situated  directly  over 
the  pan.  Due  to  the  size  of  the  enclosure,  only  a portion  of  the  curtained  area  was  modelled.  To 
reduce  the  number  of  grid  volumes  by  a further  factor  of  two,  it  was  assumed  that  the  fire  was 
symmetric  about  a plane  passing  through  the  fire  center  and  parallel  to  the  draft  curtains. 

Initial  modelling  efforts  concentrated  on  the  gas  flow  interacting  with  the  draft  curtains  and  assumed 
that  the  open  trusses  would  have  minimal  impact  on  the  gas  flow.  Heat  flow  to  the  ceiling  was 
initially  ignored  due  to  the  low  expected  temperatures  at  the  ceiling.  Later  calculations  included  a 
conducting  ceiling.  Thermocouple  measurements  prior  to  the  fire  indicated  that  there  was  only  a small 
thermal  gradient  present  in  the  hangar. 

Several  grid  volume  sizes  were  used  in  an  effort  to  determine  what  effects,  if  any,  were  caused  by  the 
grid  volume  coarseness.  A section  view  of  the  three  dimensional  grid  used  for  some  of  the  runs 
discussed  later  in  the  paper  is  displayed  in  figure  13.  The  grid  structure  is  nonuniform  in  an  attempt 
to  maximize  the  grid  volumes  in  regions  where  the  temperature  gradients  are  largest  or  where  small 
structures  are  modelled.  Details  of  the  grid  dimensions  are  given  in  table  7. 

Since  the  curtained  area  was  not  modelled  in  its  entirety,  pressure  boundary  conditions  were  used  to 
bound  the  modelling  regime  at  gas  interfaces.  At  a pressure  boundary  condition,  it  is  assumed  that 
the  inward  or  outward  flow  is  continuous  with  upstream  values  of  temperature,  scalars,  and 
turbulence  used  if  flow  is  outward  and  ambient  values  used  if  flow  is  inward.  The  pressure 
boundaries  were  kept  a minimum  of  twelve  grid  volumes  away  from  the  edge  of  the  fire  in  order  to 
minimize  any  error  introduced  by  these  boundaries. 

One  of  the  major  decisions  in  modelling  the  fire  induced  plume  was  the  choice  of  turbulence  model 
used  in  the  calculation.  For  this  work,  the  k-e  model  was  used.  This  model  has  a number  of  user 
specified  parameters  which  affect  the  ultimate  viscosity  and  heat  transfer  in  the  plume.  A number  of 
investigators  have  proposed  sets  of  parameters  for  the  k-e  model  [20,21].  The  default  values  used  in 
the  CFD  model,  HARWELL  FLOW3D,  for  Cj  , C2  and  are  the  parameters  first  proposed  by 
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Launder  and  Spalding  [17].  For  buoyant  plumes,  an  additional  term,  C3,  to  account  for  the 
production/destruction  of  turbulence  due  to  buoyancy  may  be  included.  Additional  changes  have  been 
proposed  for  some  of  the  constants  in  order  to  match  existing  plume  correlations  with  CFD 
modelling.  Several  sets  of  parameters  for  the  k-e  model,  see  table  6,  were  used  in  modelling  the 
hangar  fire. 

In  modelling  the  fire  for  these  calculations,  it  was  assumed  that  35  percent  of  the  total  heat  release 
rate  was  radiated  from  the  plume  at  the  fire  source.  The  remaining  65  percent  was  then  used  as  the 
convective  heat  release  rate  which  provided  the  energy  to  drive  the  fire  plume  flow.  The  35  percent 
radiative  loss  rate  has  been  the  standard  value  found  in  many  computer  fire  models.  For  this  fire,  as 
mentioned  earlier,  the  radiative  loss  rate  may  range  from  20  to  35  percent  [5]. 

The  plume  centerline  temperature  was  determined  in  the  following  way.  Unlike  the  thermocouple 
measurements  fixed  near  the  ceiling  which  maintained  roughly  constant  temperatures  throughout  the 
first  150  seconds  of  the  experiment,  the  thermocouples  suspended  at  the  geometrical  plume  center  did 
not  reach  constant  temperatures  until  about  100  seconds  into  the  experiment.  Visual  inspection  of  the 
videotape  suggests  that  the  swaying  of  the  thermocouple  tree  and  some  plume  asymmetry  may  be 
responsible  for  the  lack  of  constant  temperatures  early  in  the  measurement.  Therefore,  the  plume 
centerline  temperature  was  obtained  by  averaging  the  thermocouple  measurements  taken  between  100 
and  120  seconds  which  should  give  a "steady  state"  plume  centerline  temperature  profile  for  the  fully- 
involved  pan  fire. 

Figure  14  displays  a comparison  of  the  calculated  and  measured  plume  centerline  temperatures  as  a 
function  of  height  above  the  floor  for  each  set  of  turbulence  constants.  The  calculations  are  based  on 
a time  of  45  seconds  after  the  start  of  the  fire  with  time  zero  equal  to  8250  kW  (7810  Btu/s).  The 
standard  k-e  constants  with  no  buoyancy,  k34,  consistently  underpredict  the  centerline  temperature 
and  while  the  shape  of  the  temperature  profile  in  the  lower  plume  is  simulated,  the  upper  plume  is  not 
well  represented.  Including  the  buoyancy  term,  k35,  causes  the  model  to  overpredict  the  centerline 
temperatures  and  produces  a much  poorer  agreement  with  height.  Modifying  two  of  the  standard 
constants,  calculations  k36  and  k38  provide  better  agreement  with  the  temperature  near  the  ceiling, 
but  poor  agreement  with  the  overall  shape  of  the  centerline  temperature. 

A comparison  of  the  experimental  data  to  the  McCaffrey  and  Heskestad  correlations  [18,19]  for 
plume  centerline  temperature  are  given  in  figure  15.  The  correlations,  which  are  based  on 
experimental  measurements  conducted  in  smaller  spaces  and  axisymmetric  conditions  do  not 
adequately  represent  the  present  data.  The  shape  of  the  centerline  profile  does  not  provide  for  the  flat 
temperature  profile  near  the  ceiling  observed  in  the  experiment.  Calculations  k34  - k38  are  also 
compared  with  the  McCaffrey  result  in  figure  16.  Calculations  based  on  the  turbulence  constants 
used  for  k36  and  k38  provide  the  best  agreement  with  the  McCaffrey  correlation. 

One  of  the  major  difficulties  in  modelling  the  hangar  is  that  the  roof  area  is  not  only  partitioned  by 
draft  curtains  but  also  contains  a number  of  beams,  as  mentioned  earlier,  which  could  act  as 
additional  smoke  blockages  and  could  account  for  at  least  part  of  the  flat  temperature  profile  observed 
near  the  ceiling.  Figure  17  displays  a comparison  of  a CFD  calculation  using  the  turbulence  constants 
of  k36  with  the  experimental  data.  For  this  calculation,  the  grid  volumes  were  made  smaller  near  the 
ceiling  such  that  the  3.7  m (12  ft)  draft  curtain  was  described  by  12  grid  volumes  rather  than  three  as 
was  done  in  the  earlier  calculations.  The  adiabatic  assumption  was  relaxed  and  the  ceiling  was 
assumed  to  be  conducting  and  made  of  steel  which  was  assumed  to  be  of  thickness  of  0.15  m (0.5  ft). 
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While  the  inside  of  the  ceiling  was  corrugated  steel,  the  makeup  of  the  ceiling  above  the  corrugation 
was  not  obtainable.  Over  the  duration  of  the  calculation,  30  s,  the  material  above  the  steel  was 
unimportant.  It  should  be  noted  that  the  earlier  calculations  were  run  for  45  s while  this  calculation 
was  run  for  30  s.  The  shorter  run  time  was  used  for  the  larger  number  of  grid  volumes  after 
verilying,  using  the  smaller  number  of  grid  volumes,  that  the  calculation  had  essentially  reached 
steady  state  at  30  s.  Also  included  in  the  calculation  were  large  I-beams  which  ran  perpendicular  to 
the  draft  curtain  at  the  ceiling  and  effectively  converted  each  ceiling  compartment  into  a square.  The 
additional  framing  trusses  near  the  ceiling  and  the  additional  I-beams  running  perpendicular  and 
parallel  to  the  draft  curtain  at  an  elevation  about  24  m (79  ft)  above  the  floor  were  not  modelled. 
These  modifications  provided  an  improvement  to  the  calculation  with  the  predicted  temperature  closer 
to  the  measured  temperature  near  the  ceiling  and  the  deviation  in  shape  away  from  the  measured 
temperature  being  smaller  at  the  lower  elevations. 

The  observation  that  both  the  plume  correlations  and  the  CFD  computation  suggest  that  the  plume 
temperature  should  decrease  with  increasing  height  from  the  floor  rather  than  remain  a constant  as 
was  measured  suggests  that  the  experiment  or  the  building  geometry  needs  to  be  re-examined.  While 
the  draft  curtains  and  ceiling  beams  may  play  some  role  in  the  temperature  remaining  constant  near 
the  ceiling,  they  should  not  provide  much  of  an  effect  at  heights  equal  to  one-half  the  ceiling  height 
for  a building  of  this  size  at  times  early  in  the  experiment.  The  explanation  for  the  temperature 
variation  with  height  may  be  experimental  in  nature.  One  possibility  is  that  the  plume  is  not 
maintaining  a completely  vertical  rise.  In  this  case,  thermocouples  located  at  distances  between  15 
and  25  m (49  and  82  ft)  above  the  floor  could  be  sampling  parts  of  the  plume  that  are  cooler  than  at 
the  centerline  of  the  plume.  Depending  on  how  much  the  plume  has  swayed,  a constant  temperature 
profile  with  height  could  be  experimentally  produced  using  a single  string  of  thermocouples  located  at 
the  geometrical  center  of  the  fire  as  was  done  in  this  experiment.  Based  on  plume  correlations,  the 
temperature  half-width  should  be  1.9  m (6.2  ft)  at  the  15  m (49  ft)  height  above  the  floor.  Over  a 15 
m (49  ft)  length,  for  the  plume  to  sway  1.9  m (6.2  ft)  requires  an  angle  of  about  7 degrees  as 
measured  from  the  floor.  Closer  to  the  ceiling,  the  temperature  half-width  for  the  plume  would  be 
larger  and  the  averaging  effect  would  be  less  noticeable.  Unfortunately,  additional  thermocouples 
were  not  available  to  check  the  plume  symmetry  at  lower  elevations.  Visual  inspection  of  the  plume 
using  the  video  playback  suggests  that  the  plume  is  fairly  vertical  but  visual  observations  are  line  of 
sight  dependent  and  particularly  suspect  when  small  angles  are  involved.  Inspection  of  the  flame  in 
figure  5 suggests  that  there  is  some  kinking  of  the  flame  which  may  continue  upward  until  the  plume 
impinges  at  the  ceiling.  If  the  thermocouples  are  sampling  across  the  plume  diameter  rather  than  the 
plume  centerline,  lower  average  temperatures  will  be  obtained. 

A second  possibility  would  be  the  sway  of  the  thermocouple  tree  in  response  to  the  turbulence  of  the 
plume  thus  again  sampling  a lower  average  temperature  than  the  maximum  centerline  temperature. 
The  raw  data  shows  a large  variation  of  measured  temperature  with  time  indicating  the  unsteady 
nature  of  the  measurement.  Inspection  of  the  videotapes  indicate  that  the  thermocouple  tree  was 
swaying  in  response  to  the  fire  turbulence.  Whether  these  temperature  variations  were  in  part  a 
response  to  the  thermocouple  tree  swaying  in  the  plume  cannot  be  determined  using  the  limited 
instrumentation  associated  with  the  plume. 

A third  possibility  would  be  a pre-existing  temperature  stratification  which  would  slow  the  plume, 
inhibit  entrainment  and  possibly  produce  a constant  temperature  profile.  The  measured  temperature 
stratification  is  at  most  2 °C;  hence,  the  plume  should  not  have  lost  enough  buoyancy  to  produce  the 
observed  experimental  effect. 
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To  summarize  the  results  of  this  comparison,  the  CFD  calculation  using  either  k36  or  k;38  turbulence 
constants  provide  a fair  representation  of  the  plume  centerline  temperature  at  the  ceiling.  The  more 
detailed  grid  volume  CFD  calculation  provides  a better  representation  of  the  plume  centerline 
temperature  near  the  ceiling  than  either  the  earlier  CFD  calculations  or  the  two  plume  correlations. 

Temperature  measurements  along  the  ceiling  between  the  plume  center  and  the  draft  curtain  were 
fairly  constant  in  the  experiment  with  the  average  temperature  in  the  interval  between  30  - 50  seconds 
only  decreasing  by  about  4 °C  for  measurements  0.15  m (0.5  ft)  beneath  the  ceiling.  The  CFD 
calculations  at  39  seconds  using  the  more  detailed  grid  volumes  at  0.15  m (0.5  ft)  beneath  the  ceiling 
gave  a dip  of  14  °C  as  shown  in  figure  18.  At  12.6  m (41.3  ft),  which  is  beyond  the  draft  curtain, 
the  temperature  drops  substantially  for  both  the  experimental  measurement  and  the  calculation  with 
the  calculated  temperature  in  good  agreement  with  the  experiment  at  this  location. 


6.0  RESULTS  AND  CONCLUSIONS 

1.  Inside  the  plume  region,  within  the  curtained  area,  both  the  FPETOOL  and  DETACT-QS 
computer  programs  significantly  underpredict  the  gas  temperatures.  The  correlations  of  Heskestad 
and  McCaffrey,  and  the  calculation  using  the  maximum  temperature  predicted  by  LAVENT  also 
underpredicted  the  gas  temperatures  but  were  closer  to  the  measured  values. 

2.  The  CFD  calculation  using  the  modified  k-e  parameter  provided  the  best  agreement  with 
experiment. 

3.  For  large  spaces,  a hot  gas  transport  time  delay  should  be  incorporated  into  DETACT-QS, 
FPETOOL,  and  LAVENT. 

4.  Established  plume  theory  correlations  such  as  those  of  McCaffrey  for  plume  temperature  and 
velocity  in  the  r/h  < 0.2  region  need  to  be  incorporated  into  the  three  computer  fire  models  studied, 
since  for  high  bay  spaces  the  sprinkler  or  detector  of  interest  would  normally  be  located  in  this 
region. 

5.  Reasonable  agreement  was  found  between  the  field  model  calculations  using  modified  k-e 
constants  and  the  measured  centerline  plume  temperatures  near  to  the  ceiling.  Further  from  the 
ceiling,  the  agreement  between  the  calculations  and  the  measurement  was  not  as  good. 

6.  Reasonable  agreement  was  found  between  the  CFD  calculation  and  temperature  measurements 
along  the  ceiling.  The  temperature  jump  across  the  draft  curtain  can  be  qualitatively  reproduced  by 
the  CFD  calculation. 
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Table  1. 


Heat  Release  Rate  as  a Function  of  Time 

Time  (s) 

Energy  (kW) 

0 

0 

5 

4100 

8 

8250 

200 

8250 

Table  2. 


DETACT-QS  Input  Values 

Height  of  Ceiling  Above  Fuel  (m) 

30.4 

Distance  of  Detection  from  Axis  of  Fire  (m) 

0,1.5,3.0,4.6,6.1,12.5 

Initial  Room  Temperature  (°C) 

22 

Detector  Activation  Temperature  (°C) 

999 

Detector  Response  Time  Index  (RTI)  (m-s)^^^ 

33.4 

Rate  of  Heat  Release 

See  table  1 
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Table  3. 

FPETOOL  Input  Values 


ROOM 

SPACIAL  DIMENSIONS 

Ceiling  Height  (m) 

30.4 

Room  Length  (m) 

81 

Room  Width  (m) 

12.5 

LINING  MATERIALS  (CEILING  AND  WALL) 

Material 

Thermal  Conductance 

steel 

(kW/m-K) 

0.043 

Density  (kg/m“) 

Specific  Heat 

7.85E  + 03 

(kJ/kg-K) 

0.48 

Thickness  (mm) 

100 

SPRINKLER 

RTI  (m/s)*^- 

33.4 

Response  Temperature  (°C) 

999 

Radial  Distance  (m) 

SPRINKLER  IS  NOT  A SIDEWALL 

0,1.5,3.0,4.6,6.1,12.5 

FIRE* 

Fire  Height  (m) 

0.15 

Heat  of  Combustion  (kJ/kg) 

30,450 

Rate  of  Heat  Release 

see  table  1 

VENTS 

DOOR/WINDOW 

Opening  Height  (m) 

26.8 

Opening  Width  (m) 

162 

Sill  Height 

0 

*No  smoke  detectors,  HVAC,  or  program  halts  defined. 

' Fire  parameters  listed  are  those  differing  from  defaults  program. 
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Table  4. 


LAVENT  Input  Values 

ROOM  PROPERTIES 

Ceiling  Height  ( m) 

30.4 

Room  Length  ( m) 

12.5 

Room  Width  ( m) 

81 

Number  of  Vents 

0 

Curtain  Length  ( m) 

162 

Height  for  bottom 

of  Curtains  (m) 

26.8 

PHYSICAL  PROPERTIES 

Ambient  Temperature  (K) 

295 

Material 

bare  metal  deck 

Material  Thickness  (m) 

0.1 

OUTPUT  PARAMETERS 

Final  Time  (s) 

200 

Output  Interval  (s) 

10 

FUSIBLE  LINK  PROPERTIES 

Radial  Distance  (m) 

0,L5,3.0,4.6,6.1 

Distance  Below  Ceiling  (m) 

0.15,0,30,0.91,1.5 

Radial  Distance  (m) 

12.5 

Distance  Below  Ceiling  (m) 

0.025,1.2,2.4,3.7 

Radial  Distance  (m) 

0 

Distance  Below  Ceiling  (m) 

0.025 

RTl  (m-s)*^2 

33.4 

Fuse  Temperature  (K) 

999 

FIRE  PROPERTIES 

Fire  Height  (m) 

0.15 

Fire  Diameter  (m) 

3.088 

Fire  Output  as  Function  of  Time 

see  table  1 

SOLVER  PROPERTIES 

program  defaults 
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Table  5. 


Numerical  Trials 

Trial 

Cl 

C2 

C3 

C, 

Prandtl  Number 
For  Enthalpy 

CAPPA 

K 34 

1.44 

1.92 

0 

0.09 

0.9 

0.419 

K 35 

1.44 

1.92 

1 

0.09 

0.9 

0.419 

K 36 

1.44 

1.92 

1 

0.18 

0.85 

0.419 

K 38 

1.44 

1.92 

1 

0.15 

0.85 

0.419 

Table  6. 

FLOW  3D  INPUTS 
MODEL  PHYSICS 
Turbulent  Flow  (k-e  Model) 

Compressible  Flow 
Heat  Transfer 
Buoyant  Flow 
Transient  Flow 

GRID  VOLUME  DIMENSION 


21  X 28  X 90  grid  volumes 


x: 

13  * 0.6858  m 

8 * 1.3716  m 

y- 

10  * 1.3716  m 

8 * 0.6858  m 

10  * 1.3716  m 

z: 

40  * 0.1524  m 

30  * 0.6096  m 

20  * 0.3048  m 

AMBIENT  VARIABLES 


Velocity  (m/s) 

0 

Gauge 

Pressure  (Pa) 

0 

Density  (kg/m^) 

1.2 

Viscosity  (kg/m-s) 

1.91E-05 

K(J) 

l.OOE-10 

e(J/s) 

l.OOE-11 

Temperature  (K) 

295 

BOUNDARY  CONDITIONS 

Pressure  Boundary  (Sides) 

Solid  (floor  and  non-conducting  ceilings) 

Conducting  Solid  (for  conducting  ceilings) 
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Table  7. 


Plume  Symmetry 

Mean  Temp  °C  (32-200  s) 

Standard  Deviation 

Location 

Depth  Below 
Ceiling  (m) 

Test  1 

Test  2 

Test  3 

Test  1 

Test  2 

Test  3 

East  3.0  m 

0.15 

63 

59 

64 

7 

6 

7 

from  pan 

0.3 

60 

55 

61 

7 

6 

7 

centerline 

0.91 

63 

61 

65 

6 

4 

6 

1.5 

61 

65 

T : 

5 

7 

West  3.0  m 

0.15 

53 

59 

57 

4 

6 

5 

from  pan 

0.3 

50 

55 

54 

4 

6 

5 

centerline 

0.91 

52 

58 

57 

4 

5 

4 

1.5 

49 

54 

54 

4 

5 

6 
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Table  8. 


Mean  Gas  Temperatures  (32-100  s) 

Depth  Below 

T 

measured 

Standard 

Location 

Ceiling  (m) 

(°C) 

Deviation 

Pan 

Centerline 

0 m 

0.025 

64 

(6) 

East  1.5  m 

0.15 

61 

(6) 

(r/h  = 0.05) 

0.3 

59 

(6) 

0.91 

63 

(5) 

1.5 

62 

(5) 

East  3.0  m 

0.15 

59 

(6) 

(r/h  = 0.1) 

0.3 

55 

(6) 

0.91 

61 

(4) 

1.5 

61 

(5) 

East  4.6  m 

0.15 

57 

(5) 

(r/h  = 0.15) 

0.3 

55 

(7) 

0.91 

59 

(4) 

1.5 

57 

(5) 

East  6.1m 

0.15 

59 

(4) 

(r/h  = 0.20) 

0.3 

59 

(4) 

0.91 

64 

(6) 

1.5 

63 

(5) 

East  12.5  m 

0.025 

33 

(4) 

(r/h  = 0.41) 

1.2 

35 

(3) 

2.4 

42 

(7) 

3.7 

34 

(5) 
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Table  9. 


Hot  Gas  Transport  Time 

Radial  Distance 

From  Pan  Centerline  (m) 

Time  to  Thermocouple 
Response  (s) 

0 

20 

East  1.5 

20 

East  3.0 

20 

East  4.6 

20 

East  6. 1 

20 

East  6.2  (under  draft  curtain) 

20 

East  12.5 

20 

East  18.9  (under  draft  curtain) 

44 

West  3.0 

West  5.5 

20 
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Figure  1 . Picture  of  hangar 
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Figure  3.  Ceiling  of  hangar 
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f-igure  5.  Picliiro  of  fire  at  steady  state 
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Temperature  fC) 


Time  (s) 


Figure  7.  Measured  vs.  predicted  gas  temperatures  at  a radial  distance  of  1.5  m from 
the  fire  centerline,  r/h  = 0.05.  d = distance  below  the  ceiling 
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Temperature  (°C) 


Figure  8.  Measured  vs,  predicted  gas  temperatures  at  a radial  distance  of  6.1  m from 
the  fire  centerline,  r/h  = 0.20.  d = distance  below  the  ceiling 
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DRAFT  CURTAIN 


CEIUNG  JET 


Figure  9.  Illustration  of  plume  and  ceiling  jet  regions  of  a typical  fire 
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Temperature  (°C) 


Time  (s) 


Figure  10.  Temperature  distribution  within  die  plume  measured  at  radial  distances  (r)  of 
0,  1.5,  3.0,  4.6,  and  6.1  m from  the  fire  centerline,  0.15  m below  the 
ceiling.. 
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Temperature  (°C) 


Figure  11.  Measured  vs.  predicted  gas  temperatures  at  a radial  distance  of  6.1  m from 
fire  centerline.  Temperatures  predicted  using  a 20  percent  higher  heat  release 
rate  of  9900  kW 
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Temperature  (°C) 


Figure  12.  Temperature  at  4.6  m radial  distance  from  fire  centerline  with  time  shift, 
r/h  = 0.15  d = distance  below  the  ceiling 
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Figure  13.  Two  dimensional  section  of  a three  dimensional  grid  used  for  the  CFD  calculations 
representing  a portion  of  the  fire  test  curtained  area  and  the  east  adjacent  curtained 
area 
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30.5  m 


Height  Above  Floor  (m) 

Figure  14.  Comparison  of  CFD  calculations  using  four  different  sets  of  k-e  constants  to 
the  measured  plume  centerline  temperatures 
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Figure  15.  Comparison  of  the  Heskestad  and  McCaffrey  plume  correlations  to  the  measured 
plume  centerline  temperatures 
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Figure  16.  Comparison  of  calculations  using  four  different  sets  of  k-e  constants  to  the 
plume  correlation  of  McCaffrey 
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Figure  17.  Comparison  of  a CFD  calculation  using  the  turbulence  constants  of  k;36  and  a 
more  detailed  grid  with  the  measured  plume  centerline  temperatures 
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Figure  18.  Comparison  of  the  CFD  calculated  ceiling  temperatures  using  a detailed  grid 
and  the  k36  turbulence  constants  to  the  measured  temperatures  at  various 
radial  distances  from  the  fire  centerline  and  0.15  m below  the  ceiling. 
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1.0  INTRODUCTION 


NASA  hss  numerous  large,  high  bay  spaces  that  are  used  to  perform  a variety  of  different  functions. 
Many  of  these  functions  are  critical  to  NASA  operations.  Examples  of  high  bay  spaces  at  NASA 
include  those  used  for  clean  rooms,  shuttle  simulators,  assembly/storage,  vacuum  and  vibration 
chambers,  vehicle  assembly,  and/or  testing  facilities  with  payloads.  These  spaces  represent  some  of 
the  most  difficult  fire  protection  challenges  in  that  detection  of  a fire  in  a large  space  may  be  delayed 
due  to  the  distance  smoke  and  products  of  combustion  must  travel  to  reach  the  detector,  the  large 
amount  of  ambient  air  for  smoke  dilution,  the  high  dollar  value  of  these  spaces  and  the  low  damage 
threshold  of  a clean  room.  Some  of  these  spaces  also  involve  forced  airflow. 

As  was  shown  in  section  A of  this  report,  current  computer  fire  models  do  not  provide  accurate 
temperature  predictions  in  a high  bay  space  where  the  region  of  interest  is  within  the  plume.  Design 
for  fire  protection  in  these  spaces  at  present  relies  upon  "best  estimate"  guess  work.  Therefore,  it  is 
important  for  NASA  to  develop  a methodology  by  which  to  determine  fire  protection  options  for  these 
spaces.  The  methodology  must  be  in  a suitable  form  to  be  incorporated  into  the  NASA  fire  protection 
manual  and  be  easily  usable  by  the  field  engineers.  This  project  represents  the  first  step  toward 
developing  such  a methodology.  An  existing  NASA  high  bay  space  was  modelled  using  a field  model 
to  determine  how  the  existing  fire  detection  devices  would  respond  to  various  size  fires  in  the  space. 

The  CFD  model,  HARWELL  FLOW3D  [1],  used  to  model  the  NASA  high  bay  space  has  already 
been  successfully  applied  to  the  analysis  of  fire  problems  in  both  forced  and  unforced  flow  situations 
at  NIST.  Compartment  fires  with  single  and  three  room  configurations  were  successfully  modeled  [2] 
as  part  of  a validation  project,  section  A of  this  report  extended  the  validation  efforts  to  high  bay 
spaces.  The  model  was  used  to  successfully  simulate  an  experiment  [3]  used  to  determine  the  impact 
on  the  ceiling  Jet  flow  of  beams  on  a flat  ceiling.  Finally,  the  model  has  been  used  to  analyze  and 
solve  a fire  problem  associated  with  the  Navy’s  19F4  Fire  Fighting  Trainer  where  winds  blowing  in 
from  the  ocean  adversely  impacted  the  operation  of  the  trainer  [4]. 

A clean  room  at  Goddard  Space  Flight  Center  in  Greenbelt,  MD,  was  chosen  to  be  analyzed  for 
detector  placement  based  on  the  premise  that  the  forced  air  flow  through  the  room  presented  the 
possibility  of  adversely  effecting  the  detection  of  smoke  from  a fire.  Because  of  the  air  flow  in  the 
room,  models  such  as  DETACT-QS,  FPETOOL,  and  LAVENT  [5-7]  used  for  the  earlier  analysis  are 
not  appropriate  for  use  in  this  space.  These  computer  fire  models  assume  that  the  plumes  are 
axisymmetric,  neglect  flow  velocities  in  both  the  upper  and  lower  layers  and  make  the  assumption  that 
an  upper  layer  is  formed  regardless  of  fire  size  and  flow  conditions.  The  air  flow  in  this  clean  room 
may  have  sufficient  velocity  to  negate  all  of  these  assumptions  for  small  fires.  The  purpose  of  this 
analysis  will  be  to  study  the  impact  of  fire  size  on  the  distribution  of  hot  gas  and  detector  activation  in 
the  clean  room. 

2.0  NASA  GODDARD  HIGH  BAY  CLEAN  ROOM 

The  Goddard  clean  room  is  the  largest  known  laminar  flow  clean  room  in  the  world.  The  space 
measures  38.1  x 30.5  x 27.4  meters  (125  x 1(X)  x 90  ft)  in  dimension.  The  room  is  rated  as  a class 
10  clean  room  and  is  large  enough  to  accommodate  two  full  size  shuttle  payloads  simultaneously. 

Two  decks  are  located  on  the  side  of  the  room  opposite  the  air  intake  wall.  Two  35  ton  cranes  with 
hook  heights  of  21  and  24.4  m (69  and  80  ft)  are  mounted  on  beams  near  the  ceiling.  HEPA  filters 
are  used  to  remove  particulate  matter  at  the  air  intake  wall.  A perforated  return  air  grille  is  used  for 
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the  outlet  wall  which  is  located  opposite  the  intake  wall.  The  laminar  air  flow  in  the  room  is  uniform 
from  floor  to  ceiling  in  a horizontal  direction  and  is  maintained  at  a rate  of  0.508  m/s  (1.67  ft/s)  by  a 
computer  operated  mechanical  control  system.  The  concrete  floor  is  covered  by  a polyurethene 
coating.  The  clean  room  is  used  for  assembly  and  testing  of  space  flight  hardware. 

A diagram  of  the  room  showing  the  gridding  is  given  in  figure  1.  The  solid  decks  are  shown  as 
rectangles  extending  out  from  the  right  side  of  the  grid.  The  cranes  are  approximated  by  the  two 
squares  inserted  in  the  upper  right  portion  of  the  grid.  Figure  1 represents  one  vertical  plane 
extending  from  the  floor  to  the  ceiling  of  the  clean  room.  The  grid  is  actually  a three  dimensional 
grid  of  size  25  x 20  x 24.  The  large  grid  volumes  are  cubic  with  each  side  having  a length  of  1.524 
m (5  ft)  and  the  smaller  grid  volumes  have  a vertical  length  of  0.762  m (2.5  ft)  with  a square  base  of 
1.524  m (5  ft)  on  a side. 

The  smoke  detectors  on  the  outlet  wall  are  spaced  every  3.66  m (12  ft)  in  both  the  horizontal  and 
vertical  direction.  Since  the  two  decks  essentially  compartmentalize  the  outlet  wall  into  three  spaces, 
the  initial  detector  design  required  the  placement  of  at  least  one  row  of  detectors  in  each  space. 

There  are  no  detectors  on  the  ceiling. 

3.0  ANALYSIS 

3.1  FIRE  WITH  FORCED  AIRFLOW 

Smoke  flow  from  a fire  will  be  affected  by  both  the  fire  size,  which  will  impact  buoyancy,  and  the 
forced  airflow  in  the  room.  Large  fires,  which  will  create  strong  buoyant  forces,  should  produce  a 
plume  which  will  propel  most  of  the  smoke  toward  the  ceiling  and  form  a ceiling  layer.  Small  fires 
will  have  weak  buoyant  forces  which  may  result  in  very  little  upward  movement  before  the  ambient 
air  flow  causes  the  smoke  to  exit  the  room.  The  horizontal  smoke  spread  perpendicular  to  the  flow 
velocity,  may  be  reduced  with  a small  fire,  hence  detector  spacing  for  small  fires  may  be  more 
critical  than  for  large  fires. 

In  the  model,  the  fire  was  placed  three  grids  off  of  the  floor  and  in  the  far  center  of  the  room  near 
the  inlet  wall.  The  location  of  the  fire  off  the  floor  was  suggested  by  the  raised  assembly  areas  which 
were  in  active  use  when  the  clean  room  was  visited.  The  fire  was  modelled  by  uniformly  distributing 
the  heat  release  rate  in  four  adjacent  grid  volumes.  Due  to  the  geometrical  simplicity  of  the  room,  it 
was  felt  that  varying  the  location  of  the  fire  was  not  necessary.  The  default  values  for  the  k-e 
turbulence  model  were  used  in  these  calculations. 

The  model  calculations  give  the  location  of  the  heated  gas  but  do  not  explicitly  track  smoke  particles. 
A connection  between  the  calculated  hot  gas  distribution  and  the  smoke  distribution  is  therefore 
required  in  order  to  use  the  hot  gas  distribution  to  track  the  smoke  particles.  It  is  assumed  in  this 
analysis  that  the  smoke  distribution  is  coincident  with  the  distribution  of  hot  gas.  This  assumption  is 
supported  by  the  experiments  of  Heskestad  and  Delichatsios  [8]  which  showed  that  for  a fire  under  a 
flat  smooth  ceiling,  the  optical  density  (smoke  density)  is  correlated  with  the  local  rise  in  temperature. 
While  there  was  no  forced  airflow  present  in  these  experiments,  the  smoke  did  follow  the  gas  flow 
produced  by  the  buoyant  force  of  the  fire.  Unless  there  is  substantial  deposition  of  smoke  it  would  be 
expected  that  in  forced  flow  situations,  the  smoke  would  again  follow  the  gas  flow  and  that  the  smoke 
density  would  correlate  with  the  local  temperature  rise.  Therefore,  for  the  analysis  used  here,  it  was 
assumed  that  the  excess  temperature  plots  were  sufficient  to  locate  the  smoke  particles. 
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Three  fire  sizes,  32  MW,  400  kW,  and  40  kW  (30000  Btu/s,  380  Btu/s,  38  Btu/s)  were  used  to 
investigate  the  impact  of  fire  size  on  detectability.  The  32  MW  (30000  Btu/s)  fire  was  chosen  as  it 
represents  the  maximum  heat  release  rate  for  the  default  fire  in  the  model,  LAVENT.  The  other  fire 
sizes  represent  smaller,  more  likely  fires  for  the  clean  room  with  the  40  kW  (38  Btu/s)  fire 
approximating  the  size  fire  found  in  a wastepaper  basket.  Contour  plots,  in  a vertical  plane  passing 
through  the  fire,  of  excess  temperature  (temperature  above  ambient)  are  given  in  figures  2,  3,  and  4 
for  the  three  fires.  The  darker  shaded  areas  represent  higher  temperature  areas  which  should  contain 
higher  smoke  concentrations.  Clear  white  areas  will  contain  no  smoke.  The  small  values  for  the 
excess  temperatures  are  due  to  the  large  grid  volumes  used  to  model  these  very  small  fires  in  this 
large  space. 

From  the  figures,  it  can  be  seen  that  the  larger  two  fires  produce  sufficient  buoyancy  to  allow  the 
heated  gas  to  reach  the  ceiling  and  hence  may  be  detected  by  at  least  some  of  the  wall  mounted 
detectors  near  the  ceiling.  For  the  largest  of  the  fires,  there  is  potential  for  plume  driven  circulation 
patterns  to  set  up  where  some  of  the  smoke  detectors  at  lower  elevations  may  not  detect  the  fire. 

The  small  fire  did  not  possess  enough  buoyancy  for  the  fire  plume  to  rise  to  the  ceiling  by  the  time 
the  heated  gas  exited  the  room  via  the  outlet  vents.  Hence,  the  detection  of  a small  fire  in  this  space 
will  require  wall  mounted  detectors  which  are  approximately  at  the  height  of  the  fire.  The  spread  of 
the  heated  gas  in  the  horizontal  direction  at  the  base  of  the  fire  is  shown  in  figure  5.  There  appears 
to  be  adequate  horizontal  spread  for  the  smoke  to  encounter  at  least  two  smoke  detectors. 

The  computer  fire  models  used  in  section  A assume  that  the  fire  plume  rises  vertically  to  the  ceiling. 
The  forced  ventilation  effects  on  the  fire  plume  predicted  by  the  CFD  calculations  for  small  fires 
would  not  be  predicted  by  these  fire  models. 

While  the  distribution  of  the  smoke  detectors  may  be  adequate,  the  question  still  remains  concerning 
whether  the  sensitivity  of  the  detectors  are  adequate  for  small  fires.  The  wall-mounted  detectors  in 
the  clean  room  are  photoelectric  smoke  detectors  which  have  a rating  of  5.2  - 8.1  percent  obscuration 
per  meter  (1.9  - 4.8  percent  obscuration  per  foot).  This  rating  is  related  to  the  optical  density  per 
meter  of  the  smoke  by  [9] 


Where  is  the  observed  intensity,  is  the  incident  intensity,  D is  the  optical  density  per  meter,  and 
L is  the  path  length.  D is  related  to  the  fire  by  the  relation 

D=DJ,LmlV) 


Where  is  the  mass  optical  density.  Am  is  the  mass  consumed  in  the  fire,  and  is  the  smoke 
volume. 

The  above  two  equations  must  be  evaluated  for  a given  fire  in  order  to  estimate  the  optical  density  at 
the  detectors.  The  values  of  D^,  and  Am  depend  on  the  fuel  and  also  on  the  heat  release  rate.  Fuels 
typically  have  heats  of  combustion  between  16  and  40  kJ/g  and  produce  mass  optical  densities 
between  0.12  and  1.4  m“/g  (590  and  6800  ft^/lb).  Fuels  found  in  the  clean  room  may  consist  of 
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plastics,  mylars,  metal,  and  wood.  For  the  calculations  given  below,  it  will  be  assumed  that  the 
material  heat  of  combustion  will  be  30  kJ/g  (28000  Btu/lb)  and  the  mass  optical  density  produced  will 
be  0.5  m^/g  (2400  ft^/lb).  This  will  give  a middle  range  result  with  some  materials  producing  more 
smoke  and  others  less  smoke.  The  smoke  mass  per  volume  Am/V^  will  be  estimated  using  the 
following  relation 


Am  _ Q 


where  H^.  is  the  heat  of  combustion,  Q is  the  heat  release  rate,  A is  the  area  occupied  by  the  smoke 
perpendicular  to  the  forced  airflow,  and  v is  the  velocity  of  the  forced  airflow  in  the  clean  room.  The 
area  A at  the  detector  wall  is  determined  from  the  field  model  calculations  by  assuming  that  the 
smoke  is  essentially  contained  in  the  region  at  the  wall  where  the  excess  temperature  is  greater  than 
1/3  of  the  maximum  excess  temperature. 

The  smoke  from  the  40  kW  (38  Btu/s)  fire  impacts  an  area  on  the  outlet  wall  of  approximately  140 
m-  (1500  ft^),  approximately  2/3  of  the  area  under  the  lower  deck.  Using  the  above  mass  optical 
density  and  heat  release  requires  a detector  with  a rating  of  2.3  percent  obscuration  per  meter  (0.7 
percent  per  ft)  which  is  significantly  below  the  sensitivity  of  the  existing  detectors.  The  400  kW  (380 
Btu/s)  and  32  MW  (30000  Btu/s)  fires  had  perpendicular  areas  which  approximately  equaled  the  clean 
room  wall,  840  m“  (9000  ft^).  The  400  kW  (380  Btu/s)  fire  would  require  a rating  of  3.9  percent  per 
meter  (1.2  percent  per  ft)  or  75  percent  of  the  existing  detector  rating  while  the  32  MW  fire  required 
a rating  of  2(X)  percent  obscuration  per  meter  (61  percent  per  ft). 

Both  the  40  kW  (38  Btu/s)  and  the  400  kW  (380  Btu/s)  fires  generated  estimated  smoke 
concentrations  below  the  detection  limit  of  the  existing  smoke  detectors  in  the  clean  room.  This 
conclusion  depends  on  the  smoke  following  the  heated  gas  and  being  well  mixed. 

3.2  FIRE  WITHOUT  FORCED  AIRFLOW 

The  4(X)  kW  (380  Btu/s)  fire  was  used  to  demonstrate  how  the  analysis  can  change  when  the  forced 
airflow  is  absent  in  the  clean  room.  The  heat  source  representing  the  fire  was  positioned  as  in  the 
earlier  analysis  but  a smaller  grid  was  used  to  define  the  plume  more  accurately.  Two  openings, 
simulating  open  doors,  centered  at  the  floor  on  opposite  walls  are  used  to  vent  the  room.  Figures  6 
and  7 present  the  calculated  temperature  distribution  at  times  of  50  and  100  s after  the  fire  ignition. 

As  can  be  seen  from  the  figures,  a vertical  plume  rises  from  the  fire  source,  is  deflected  by  the 
ceiling  and  forms  a ceiling  jet  along  the  ceiling.  Since  the  fire  is  not  symmetric  in  the  room,  the 
ceiling  Jet  strikes  the  near  wall  first  and  sets  up  a circulation  pattern  there.  The  smoke  is 
concentrated  in  the  ceiling  jet  and  does  not  diffuse  as  it  does  in  the  forced  airflow  environment.  The 
gas  temperature  without  forced  airflow,  see  figures  3 and  6,  is  substantially  higher  than  is  calculated 
in  the  forced  airflow  simulation.  In  these  situations,  smoke  detectors  located  at  substantial  distances 
beneath  the  ceiling  would  be  of  little  use  in  the  early  detection  of  the  fire.  Because  the  smoke  is 
concentrated  in  the  ceiling  jet,  the  requirements  on  smoke  detector  sensitivity  near  the  ceiling  would 
be  less  severe  than  in  the  forced  airflow  environment. 


52 


4.0  RESULTS  AND  CONCLUSIONS 


1.  Computer  fire  models  examined  in  section  A of  this  report  are  inappropriate  to  use  for  detector 
analysis  in  forced  air  environments. 

2.  The  design  fire  size  is  of  extreme  importance  in  structures  where  there  is  a forced  air  flow.  For 
the  NASA  Goddard  clean  room,  small  fire  plumes  would  not  reach  the  ceiling,  hence  the  distribution 
and  sensitivity  of  side  wall  detectors  becomes  important. 

3.  Large  fires  in  structures  containing  forced  air  flow  may  provide  sufficient  buoyancy  for  the  fire 
plume  to  reach  the  ceiling.  In  these  situations,  wail  mounted  detectors  may  be  shielded  by 
obstructions  from  the  smoke  and  not  detect  the  fire. 

4.  Detector  sensitivity  in  the  NASA  Goddard  clean  room  should  be  increased  for  small  fire  detection. 
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Figure  1. 


Vertical  slice  of  a typical  grid  used  for  the  NASA  clean  room  calculation 
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Figure  2.  Temperature  contours  on  a vertical  plane  through  the  fire  for  the  32  MW  fire 
calculation 
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Figure  3.  Temperature  contours  on  a vertical  plane  through  the  fire  for  the  400  kW  fire 
calculation 
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Figure  4.  Temperature  contours  on  a vertical  plane  through  the  fire  for  the  40  kW  fire 
calculation 
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Figure  5.  Temperature  contours  on  a horizontal  plane  centered  on  the  grid  directly  above  the  40 
kW  fire 
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Figure  6.  Temperature  contours  on  a vertical  plane  through  the  fire  at  50  s after  ignition  for  a 
400  kW  fire 
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Figure  7.  Temperature  contours  on  a vertical  plane  through  the  fire  at  100  s after  ignition  for  a 
400  kW  fire 
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